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THE REVOLUTION OF A DARK PARTICLE ABOUT A LUMINOUS 

CENTER 

By Edwin Bidavell Wilson 

1. Introduction. The simple problem of two bodies moving under 
their mutual gravitational attraction has doubtless long been considered as 
completely solved. Yet ever since Maxwell's formulation of his theory of 
electromagnetic action as a motion or stress of some kind in the ether, there 
has been a remote possibility that the time might come for revising the dynam- 
ical equations of motion so as to take into account the action of the ether on 
the moving masses. In fact, recent investigators of the electronic theories of 
matter have paid a great deal of attention to the interaction of matter and 
ether. To adapt their work to the problem of the revolution of bodies about 
a fixed luminous center, which is approximately the case that obtains in the 
solar system, is by no means easy. For in the first place there are diflTerent 
theories of electrons and in the second place these theories lead one to try to 
replace the older and more familiar dynamical equations by newer and as 
yet not fully settled electrodynamical equations. Thus Abraham* discusses 
in considerable detail the actions of the ether on a moving body ; but his 
equations are those of the electromagnetic field and depend to some extent 
on his particular theory. 

Recently, however, Poynting, with that directness of insight and sharp 
mechanical instinct which characterize the English investigators of physical 
science, has taken up the problem of the effect of radiation on moving masses 
and has treated it from the point of view of the transmission of momentum by 
the ether. His results, although perhaps essentially electromagnetic in that they 
are derived from the consideration of waves in the ether, are nevertheless 
well adapted to direct mechanical formulation, t The effects which he describes 

* Theorie der Eleklrizitat, vol. 2. 

t We do not wish to be interpreted as implying either that the pressure of radiation is 
necessarily connected with any electromagnetic theory or that Poynting was the first to con- 
sider the matter from otiier points of view. On the contrary the question is rather one of 
mechanics and thermo dynamics and had previously been treated as such by several authors. 
Historical references may be found in an article on " The pressure due to radiation " by Nichols 
and Hull, Astropkysical Journal, vol. 17, June, 1!)03 or The Smithsonian Seport for 1903, pp. 
1 15-138. l^oynting merely reverts to this simpler treatment and applies It anew. 

(135) 
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are three :* the direct pressure of light, the Doppler reception effect, and 
the Doppler emission effect. We will state what is meant by each. 

If a body which is under the rays of the sun be considered as dark t or 
essentially black, to use the technical term, the radiation from the sun is 
stopped and absorbed by the body. The result is an outward pressure due to 
the abstraction of momentum from the ether. As the electromagnetic radiation 
varies in its intensity according to the same law as the gravitational attraction, 
it produces a pressure which varies inversely as the square of the distance and 
hence introduces no change in the dynamical equations except an apparently 
diminished value for the gravitational action. It might seem at first sight 
that in estimating the pressure due to radiation the distance should be meas- 
ured from the radiating surface : but a moment's consideration shows that a 
spherical body like the sun could be replaced by a concentric body radiating 
with a properly altered intensity, and hence that the distance should be meas- 
ured from the center. 

The two Doppler effects are much smaller than the direct pressure of light. 
To consider the reception effect suppose that the dark body were moving 
toward the sun or any other luminous center with a velocity v. Then a 
greater amount of momentum would be absorbed per unit time than if the body 
were at rest. In particular the amount would be greater in the ratio vj V, where 
V denotes the velocity of light. In case the body falls to the sun from an in- 
finite distance this ratio is about 2/1000, which may be considered as its maxi- 
mum value in the solar system ; generally it is much smaller. If the body 
retreats, the pressure is lessened in the same ratio. Finally there is the emis- 
sion effect which arises from the fact that the dark body radiates the energy it 
receives. This force, which is "frlctional" in that it always opposes the 

♦ See his address to the Royal Society, Nature, vol. 76, pp. 90-93, November 22, 1906. 
Earlier than tliis, Poyntiug had printed an extensive article on " Radiation in the solar system : 
its eflfect on temperature and its pressure on small bodies," Transactions of the Boyal Society, 
ser. A, voi. 202 (1903), pp. 525-652 In this article there occur some minor errors which 
the author corrects in a note at t!ie end. Moreover the two Doppler effects are by no means 
clearly distinguished one from another. Finally in this earlier treatment, loc. cit., p. 649, 
the author comes upon a differential equation that he states is probably not integral)le 
This equation corresponds to our equation (5). Tlie integration which we give in terms 
of Bessel functions not only furnishes a desirable complete determination of the problem, but 
enables us to see certain other general characteristics of the motion than are obtainable from 
a less exact integration. 

t This is approximately true for the small particles in the solar system. As far as we 
know them, their albedo is very low. 
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motion, is stated by Poynting to be one-third of vj V of the direct pressure of 
light.* This does not make the emission effect merely one-third as great as 
the reception eflfect : for the velocity which enters into the determination of 
the reception eflfect is the component velocity of the motion resolved along the 
line joining the particle and the sun, whereas in the emission effect it is the 
total velocity in path that enters into the formula. 

2. The differential equations of motion and their integrals. 

These statements of Poynting's may be incorporated in the dynamical equations 
in a simple manner if polar coordinates with the sun as center be assumed for 
reference. Let — Cjr^ and A'/?"* respectively be the accelerations due to the 
gravitational attraction and the direct pressure of radiation. There will be the 
accelerations 

r'' V di Ir'' V di 

along the radius vector due respectively to the reception and emission effects. 
Further there will be the single acceleration 

_\K\_rde 
3 r'' V dt 

perpendicular to the radius. The equations of motion then become t 

O'r rdeV 0-K i K dr 



^ ' d^ ^ \dtj ~ r* 3 Vr'' dt ' 

^^ r dtV dtj' 3 Fr* dt 



The simplicity of these equations is largely due to the fortunate linear 

* As a matter of fact Poynting states, as Is well known, that higher powers of v/V 
enter into the expression ; but neglecting them does not introduce, on the assumptions 
made above, an error of more than one-flfth of one per cent in a quantity already small; 
and generally the error is much less than this. To introduce the exact expression for the 
emission effect (which is evaluated on assumptions that are not perfectly well established) 
would unnecessarily complicate the equations of motion. 

t Provided we disregard, as we may, any changes in the mass of the particle which 
might arise from its motion through the electromagnetic field. On the electronic theories 
there is such an effect, but it remains very small until the velocity of the moving body is 
within a few percent of the velocity of light. 
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dependence of the Doppler effects on the velocity of the particle. Equation (2) 
is immediately integrable into 

To find the orbit of the particle we may eliminate the time between this equa- 
tion and (1). Meanwhile in order to reduce the resulting equation to a 
simpler form, it will be well to introduce the change of variables 

(4) u=~, ^-e = cj>. 

The differential equation of the orbit then takes the particularly simple form 
d^_3dM _ 9V^ C -K _L 

Two things may be noted at this point. First, if K be allowed to ap- 
proach zero, «/) becomes infinite in such a manner that K<^ remains finite and 
d<f> = — dd. Hence (5) approaches, as it should, the form 

^ + M = const., 

ordinarily obtained in the theory of gravitation for the revolution of a planet. 
Secondly it is noteworthy that (5) is, except for the presence of the second 
member, an equation of the Bessel type. 

To bring the differential equation of the orbit into shape for easier com- 
parison with the treatises on Bessel functions the change of variable u = z^* 
will be made. Then (5) takes the final form : 

dh 1 dz /, 2»\ L 

(^> w^4>d4>^\}-j^r^-^^- 

If the second member were zero, the general solution would be 

z = AJ^{4,) +BY^(<I>) 

in terms of Bessel functions e/and Neumann functions Y.* It remains to find 
a particular solution of the complete equation to add to this. These particular 

* The notations throughout this article will conform to those of the Treatise on Bessel 
Functions by Gray and Mathews (1895). 



1907] HEVOLUTION OF DARK PARTICLE ABOUT LUMINOUS CENTER 139 

solutions, rarely given in the books, may readily be obtained by the method 
of variation of constants. In this case such a solution is 

Hence the equation of the orbit may be written as 

where the constants A and JB may be chosen so that the lower limit of inte- 
gration may take any desired value. 

As the Doppler effects are very small, the motion of the particle will be 
nearly elliptic, parabolic, or hyperbolic. To bring out more clearly the de- 
partures from these forms we may expand the solution (7) by means of the 
asymptotic expressions for the Bessel and Neumann functions. These series 
will semi-converge with extreme rapidity because of the large value of the 
argument <f) (we shall have no occasion to consider cases in which <}> is mat- 
erially less than 10* as may be seen by reference to p. 148). For such cases as 
these the function Y is inconvenient : it contains an additive multiple of J. 
Let us write * 

Y„{4>) ^ (log 2 - y) J„{4>) + Z,. (<^). 

Then the function Z satisfies the fundamental relations which t/and l^satisfy 
and which are the reduction formulasf especially useful in the manipulation 
of these functions. We have then the asymptotic formulas 

/ir<^\V2 f(2«+l)7r ■] „ . f(2n+l)7r ^^ 
\-Y) '^»("^) ^ ^ ^^4^ 4 - <^} + « ^^^\- r^ - <^|. 

/2<^\V2 f(2w+l)7r I f(2M+l)7r 1 

WJ "^"^^ " ^ 1^ 4 </>|-Psm|^ __Z__^l, 



where 



P = 1 — ^ ^ ^ ' + higher even powers ot j-, 



4,t2 _ 1 (4«2 - 1) (4w2 _ 9) (4ft2 _ 25) , . ^ , ^ ,1 

Q = -^^ - ^ ' 3!(8f)3 + '"gl^^'- °^^ P«^«^^ «f ^- 

* See Gray and Mathews, p. 40, formula 95. 

t These are: Z'„ = (n/0)Z„ — Z„ + i =Z„_i— (n/*)Z,., Z„ + i— (2n/0)Z,. +Z„_i = 0, 
<f,t.->r\ Zn — JnZn + 1 = 0, and {<t>"Zny = <f>"Zn-i. See Gray and Mathews, p. 14 and p. 16. 
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Or if we set (2n + 1) 7r/4 — <f> — ^fr„ we may write 

(-y) -4(0) = cos ^„ + -2^ sin ^„ - ^^ 2ll2<t>y cos Vr„ + . . . , 

•'* (^^)"'Z.(« = - .„ ^„ . "-Ij^ e« ^. . '"'-*»;;.-» sin t. . ■ • ■ . 

with the relations cos -</r,j_j = sin i|r„, sin i/r„_i = — cos i/r„. The terms 
omitted from (8) are of the order 10"^^ when (f> is of the order 10^. 

In terms of the functions J and Z the equation of the orbit becomes 

where C and D are constants of integration slightly different from A and Ji. 
This form of the equation is not yet in the shape most convenient for computation. 
For the integrals in the bracket change their value by a hundred percent when 
^ runs through an interval of 2'jr. It may be noted, however, that the maxi- 
mum value of the integrals, when extended over a small interval in (ft, is very 
small owing to the fact that the denominator remains of the order 10*^ and the 
numerator of order 10~^. An integration by parts will enable us to ac- 
centuate this. We have 



Hence 



/ 



•W)., = ^ + «-iffi + 6.«^ + 



with a similar result for the integral of Z2(<^)/<^^. Combining, we have for 
the bracket 

^e(<^)/^ # - M<f>)J^^ ^*'' = ^4 - ^fe + higher even powers of i* 
There appears to be no necessity of keeping terms of the order 10 ~**. Hence 

(10) u = <^*| CM<f>) + DZ,(4>) + i(ij - ^) I 
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may be taken as the approximate equation of the orbit for future compu- 
tations. For the derivative we have, from (10), the expression 

(11) «' = <^^|aJi(<^) + DZ,{4>) -L^^, 

where the small quantity 48 Lj^'' has been neglected within the brace 

If the initial values Mq, m'o, and ^o oiu, u', and <f> be substituted in (10) 
and (11), the constants C and D may be obtained. 

= l(ii!Zi(<^o)-^Z2(<^o)), 

= i(-i?Ji(<#.o) + A9J,(<^o))- 
In equations (S) let yjr„ = yjr^ = yjr. Then we have, if (2n + l)7r/4 — <^o = •^o» 

-f } '^i('^«)=«i» to - ^cos V^o - 2^p^, sin V^o + • • • . 

(^— j ^(.^o) = costo + 4sinV.o- J^(2^.cost«+ • • • • 
Hence 
Zii<t>o)Ji(.<l>) = \X-l) I *^os -f J cos -f + cos -f sin i/r ^ + sin ^r, cos l/r -i- 

+ sinV^o sinf t_^ - costocost^la^ - <^°« V^« ««« V^ 2!(i^»] • 
and 

<A(<^o) Zi(((>) = r T-^ j I - sin -fo sin -f + sin yj^Q ^os i|f ^ + cos i/ro sin 1^ ^ 
- cos to cost ^ + ^i^V^osinta^g^, + sin to sin ^ ^^^^J • 
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Hence 

BZ,(<t>o)M<f>) - IiJi{cf>o)Z,{<t>) = B (±^'''^cos(ylr, - t) 

In a similar manner 

- SZ^{<(>o)M<t>) + SJ^(<f>o)Z,(cl>) 

(I \i/2r /15 15 \ 

^) [8in(to - V^) + cos (V^o - t) (s^ - 8^j 

• / . . N / 225 105 105 \-l 

+ sin (to - ^1^) (^gj^ - 128^. - 128^^; J " 

Now let A^o — ^ = ?• Then cp — (f>Q= ^ and <^ = <^o (1 + l/^o) • Hence, by 
expanding into series for two terms and collecting, we have 

M= (P) jiZcosf + /Ssinf -i^sinf ^X + ^f sinf ^, - /S'^cosf g^, 

,. . . 240 1 r^,/l 12\ 

It should next be noted that the coefficients M and S consist of the sum of 
terms of different orders of magnitude. In fact Mq — Zi/<l>o may be considered 
as of the second order whereas 12 L/(f>l is of the fourth. So also in /S we have 
u'o of the second order (unless the orbit differs only infinitesimally from a 
circle, u' and u are of the same order) whereas 2i/0^ is of the third. If we 
rearrange the expression for u on this basis and at the same time expand 
(<^/</)o)*''* we get as the final form. 

1 ri5 , . . 24i . ^_ 15 _ . 3/ »Ai. • fc 

+ W^' + K^isinf + («o - I) cos^)r]. 
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and hence, on differentiating, 



m' = Uj cos ^ - f Mo - 4 j sin H 

1 rex V. 3 . . 69 2/ . . 21 , . . J. 

where the first, second, and third brackets denote respectively the first, second, 
and third approximations. The integration of the problem, both in exact 
and in approximate form, may therefore be considered as accomplished. 

3. Discussion of the characteristics of the motion. In the first 
place let it be assumed that at the start the areal velocity is positive. Equation 
(3) shows that as the motion goes on and increases from its initial value 0, 
the areal velocity will steadily decrease and hence by equation (4) the argu- 
ment <f> will decrease. If we allow </> to approach zero in the complete in- 
tegral (7), it appears that u is becoming infinite : for the terms 

<^V,(<^), <^*F,(<^), 4'''J,(4>)J ^, 

approach finite limits and the integral* 



^^W# 



r- 



d<t> 



diverges for <^ = 0. Hence the particle will reach the attracting center when 
^ = 0. It should be noted, however, that in an actual case the motion will be 
stopped before the particle reaches the center by the impact of the particle on 
the surface of the attracting mass. For another reason it is also valueless 
to investigate what happens to the orbit for very small values of 4>, that is, 
when the particle is near the attracting center. We have assumed in setting 
up our equations that the velocity of the particle did not become very great. 

*This follows immediately on substituting the principal values of the functions ob- 
tainable from the developments in series. See Gray and Mathews, p. 11 and p. 14. 
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Hence the particle should not be allowed to approach the center too closely. 
Hence 3 J'lijK is an upper limit for 6, and we may say that the number of rev- 
olutions that the particle can make is less than 

3 FA 



(14) jsr = 



27rK' 



This number for any given particle depends solely on the initial areal velocity 
and not on the velocity of projection or the initial distance from the center, 
except in so far as they conjointly determine the areal velocity. 

The approximate solution (12) will serve to discuss the motion for a few 
dozen revolutions from the initial position ^ = <^^. If it be assumed that the 
polar axis is so chosen that mJ = 0, that is, so that the particle starts from an 
apse, then the first approximation may be put in the form 

L 

T ^^ ~ ^ 

u = :i-^ (I + ecos ^), e= =^" • 

This is a conic section as might have been expected. In particular it is that 
conic section which the particle would described under the action of the gravi- 
tational attraction diminished by the action of the direct pressure of light and 
free from the Doppler effects. This conic I shall call the osculating conic 
of the orbit. 

To find the amount by which the particle falls in from the apse toward 
the center during the first complete revolution it is merely necessary to sub- 
tract the values of the radius vector as given by CI 2) for ^ = and ^ = — 27r. 
Hence to the second approximation we have 



^^» = -'»^rT^(^-|0' 



This quantity A rj is always negative, showing an approach to the center in all 
cases,* and varies between the limits r(,/4iVand 2rj/iV"as the osculating conic 
varies from a parabola to a circle. It is to be noticed that, if the eccentricity 
is small, the particle falls in toward the center by an amount about twice as 

• We have assumed tacitly that Z. > 0, that is, that the attraction is greater than the 
pressure of light. 



1907] REVOLUTION OF DARK PARTICLE ABOUT LUMINOUS CENTER 145 

great as that required to carry it to the center in iV revolutions if it moved in 
by the same amount during each revolution. This need not be perplexing, 
for according to the formula the particle is always falling in by a certain per- 
centage of the distance remaining (not of the original total distance) . It there- 
fore appears necessary that the particle should fall in by an amount which at 
any rate for orbits of small eccentricity is decreasing from revolution to revo- 
lution. In fact if the increment of the radius vector be found to another ap- 
proximation, it appears that such is the case.* 

As the eccentricity of the osculating ellipse seems to be of importance, it 
appears worth while to calculate the change of eccentricity during a revolution . 
Formula (13) shows that if the particle starts at an apse, it will be at an apse 
after a complete revolution provided we disregard the effect of the last ap- 
proximation. This fact will materially simplify the discussion of the eccen- 
tricity. It enables us to write the equation of the osculating ellipse after one 
revolution as 

){ = ( Wo + At(o - -ri 5— X2) cos I + y-r g— 

On calculating the eccentricity of this ellipse by the usual formula we find 

Hence it appears that the ellipse is always losing eccentricity. To speak 
roughly one may say that the orbit is losing its eccentricity from about two to 
about fifteen times as fast as it loses its perihelion distance according as the 
orbit is approximately circular or approximately parabolic. 

Equation (13) shows that although there is not much change in the angu- 
lar position of the apse, there is some. Let us assume that the apse advances 
by the amount e, so that if | = is a root of u', then ^ — — iir — e is the 
next but one. It has been seen that e is small of the order 1/iV*. The second 
term of the first approximation gives a term in e. All the terms of the second 

* If I understand Poynting aright he states In Nature, loc. cit. p. 92, second column, 
that the particle will fall in farther during the succeeding revolutions than during the pre- 
ceding ones. As Poynting very seldom makes a mistake, I have given the matter a great 
deal of attention but fail to And any error In my work. Moreover Poynting, in his earlier 
paper, obtains a formula which is interpretable in much the same way as mine. Unfortu- 
nately I cannot find that he has printed any mathematical treatment of the subject since 
then. 
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approximation are too small to contribute anything. The first and fifth terms 
of the third approximation contribute quantities independent of e. The equa^ 
tion for e therefore is 



Hence 



^-n-m e 



K*'^ + r)=^re-C^ + ^)^'^''''' 



and there is a slight advance of the apse in the direction of motion. If the 
eccentricity is very small, that is, small of an order approaching XjN^ the apse 
may advance over a considerable angle during a revolution.* 

So much for the discussion of what happens in the vicinity of the start. 
We are now in a position to take up more general considerations. From the 
behavior of the apse investigated above, it should appear that as <^ decreases 
from its initial value ^o> the apses draw further and further apart ; and con- 
versely if we follow the motion backward allowing <f) to increase, the apses will 
follow each other at intervals diminishing toward the limit ir. Now this is 
Avhat we should obtain on general considerations applied to formula (11). 
For the apses correspond to roots of u' , that is, to roots of the expression 

a T 

GJM) + ^^x(<^) - -^' 

Here C and D are definite constants. To be sure they are very small in the 
eases we have been considering, in fact of the order <^7,''-. Nevertheless for 
values of ^ sufficiently large, the third term, — 2L/(f>^, may be neglected rela- 
tively to the first two. Then the roots of u' may be investigated by the usual 
methods employed in discussing the roots of Bessel functions, f and we shall 
thereby obtain precisely the result concerning the advance of the apse that we 
have already obtained from our computations. As the value of u' becomes in- 
finite when ^ = 0, it is clear that u' must have a definite smallest root corre- 
sponding to the last apse which the orbit possesses, after which it coils in 
around the attracting center with an essentially spiral form. It is not impos- 

* It is, however, by no means certain tliat tlie orbit would Iseep Its approximately con- 
ical nature for sucli small values of e : it might already have become spiral. See infra. 
t Gray and Mathews, chap. 5. 
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sible to get a very rough estimate of where this last root occurs. For there 
can be no roots for values of 4> so small that the third term above is larger 
than the other two. If we compare merely the magnitudes of <f), (f)^, G, and D, 
we find that the condition is fulfilled when <f> = ^y^. To get any more ac- 
curate estimation seems not to be worth while. 

Similar considerations may be applied to formula (10) to discuss the roots 
of u. If the test suggested above be used, it appears that the roots of u are 
likely to be found for values of <f> of the same order of magnitude as the initial 
value 00. That this is correct follows from the approximate formula for u. 
Here, if the eccentricity e is greater then unity, there are roots of u in the 
immediate neighborhood of <f>o. If e is considerably less than unity, the 
motion of the particle may be followed back through a very large number of 
revolutions before the eccentricity shall have increased to unity. Nevertheless 
the rate of increase of e as given by (15) is tolerably fast, and it would be 
safe* to say that the first root of u would appear before <f> had increased to 
many times ^(,. Now whenever u has a root, the particle has retreated to in- 
finity. Hence it follows that not only is the number of revolutions which a 
particle can make inward toward the center a finite number JV, but the total 
number of revolutions which a particle can make in infinite time about an at- 
tracting luminous center is finite. This result is so different from that of the 
ordinary theory of attraction as to be very striking. In fact, it could hardly 
be credited if the analysis had not assured its correctness. For values of <^ 
greater than the first root of u the physical interpretation is simply that the 
particle is describing a hyperbolic type of path. As <l> gets larger and larger 
the roots of u become sepaiuted by intervals equal more and more nearly to 
TT. In other words the path is becoming more and more nearly straight. The 
analytical deductions, therefore, are seen to fulfil in every respect the con- 
ditions anticipated by the physical intuition. 

• In fact if we express C and D in terms of e„ and if we assume that the particle 
starts at an apse where <p„ has sach a value that 22(<^„) = 0, then Ji{<t>^) = and the equa- 
tion for roots of u may be written roughly in the form 

1/2 



^^|-..«2:(.o)^.W=|, ^•(*o)=(4) 



Hence we may look for roots of u in the neighborhood of 

"ft 
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In closing we may take the time to obtain numerically some of the 
quantities which we have been discussing. The case of particular interest is 
that of a particle revolving about the sun. The following data may be assumed : 

velocity of light, V=B X lO^' 
distance of sun from earth = 493 X V 
acceleration of earth in orbit = 59 X 10"* 
pressure of sunlight at earth's distance = 59 X 10'^. 

All these quantities are in the C. G. S. system. The constant C in equation 
(1) has the value 

C= 59 X 10-2 X 493-2 V^= 1-29 X lO^s. 

The constant Ji'may be determined in terms of the radius a and the density p 
of the particle as 

K: e = 59 X 10-6 ^^j2 . 4.n-p(,,359 ^ 10"^, K= fp-VrUO-^C. 

If the particle has a density about equal to that of the earth and a radius of 
about 0-Olcm, then K = 1-76 X 10*^. It should be noted that ^is increased 
by a decrease in the radius or density, and that the value of Ii here obtained 
is certainly small for such particles as meteoric dust finely divided. Suppose 
that the initial double areal velocity is 45 X 10^** which is roughly that of the 
earth. Then if a particle be started at an apse at the earth's distance from the 
sun, 

<^o = 2-3 X 10', JSr= 3-66 X 10«, e ^ 0-052. 

Hence the distance of fall from an apse will be 80- 6km during the first 
revolution. The change of eccentricity will be about one millionth of its value. 
The apse will advance about one half of one millionth of a second of arc. It 
is very doubtful whether the particle could make so many as 10' revolutions 
altogether. If the particle had a radius of one tenth the amount above as- 
sumed, the fall would be ten times as great and so would the diminution of the 
eccentricity ; the rotation of the apse, however, would be one hundred times 
as great as it was before. 

Yale University, 

New Haven, Conn. 



